Short-length storage of intense optical pulses in solid by adiabatic passage 
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We propose a novel scheme of storage of intense pulses which allows a significant reduction of 
the storage length with respect to standard schemes. This scheme is particularly adapted to store 
optical information in media with fast relaxations. 

PACS numbers: 42.50.Gy, 42.50.Md 



I. INTRODUCTION 



The propagation of two pulses in resonant interaction 
in a A- atomic medium has been widely studied during 
the last decade, in particular with application to quan- 
tum information processing (see, e.g., reviews The 
possibility of quantum information storage in the gas 
phase has been demonstrated in cold atoms Q and in hot 
ensembles Q. Achieving information storage in solid- 
state systems, which are more attractive due to their 
higher density^compactness, and absence of diffusion, has 
been pursued [J] . The main drawbacks preventing the ef- 
ficiency of storage for such solid-state materials are huge 
inhomogeneous broadenings and high rates of decoher- 
ence. For instance, in crystalline films doped with rare- 
earth metals the rates of transverse relaxations amount 
to tens of GHz 5]. In practice an efficient storage of in- 
formation requires rather large optical length such that 
even a weak loss rate will ruin it [6|]. In order to re- 
duce the inhomogeneous broadening, it was proposed in 
a number of works to use the so-called hole burning tech- 
nique 0] • It has been shown that media prepared in such 
a way allow efficient coherent population transfer in A 
systems Q. However, the hole burning technique faces a 
so far unsolved problem that leads to the reduction of the 
optical length of the samples, which will be detrimental 
in general for the efficiency of the storage. 

In the present work we show that, for the storage of 
optical information (i.e. of classical fields), the recording 
length can be dramatically reduced with the use of in- 
tense pulses. The possibility to store and retrieve optical 
information in resonant media has mainly been studied in 
the "linear approximation" with respect to the so-called 
probe field that has to be stored, i.e. for a weak probe 
pulse. It was usually assumed that the control pulse 
propagates in a medium without pulse shape change. Nu- 
merical simulation of this problem without restrictions 
on the probe intensity was performed in [9(. Analytical 
studies of the problem taking into account the group ve- 



locities of both pulses were performed in [To| in the limit 
of pulses of duration much shorter than all the relaxation 
times, but sufficiently long to allow adiabatic evolution 
during the interaction. It has been found that the length 
of information storage depends remarkably on the ratio 
between the oscillator strengths of the adjacent transi- 
tions that determine the group velocities of the pulses in 
a medium. The present work aims at a complete ana- 
lytical study of information storage in case of arbitrary 
relaxation times and intensities of the probe and con- 
trol fields. We show analytically and numerically that 
it is possible, for a proper choice of solid-state medium, 
to dramatically reduce the optical length needed to store 
intense short pulses. This would in particular make more 
efficient the use of the hole burning technique. 

The paper is organized as follows: We first describe 
the model, and analyze the propagation in the adiabatic 
limit. We next derive the conditions of storage for short 
lengths of the medium, before concluding. 
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FIG. 1: Level scheme diagram for a A-system interacting with 
two laser pulses. 



II. THE MODEL 

We consider a medium of A three-level atoms with two 
ground states 2) and a metastable excited state |3), 
of energy fiiOi, i — 1,2,3, and two laser pulses referred to 
as the probe and control pulses, coupling the transition 
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respectively 1 <-> 3 and 2^3 (see Fig. 1). The inter- 
action Hamiltonian in the rotating wave approximation 
can be written in the basis {|1), |2), |3)} as 



H 



-CI* 



2A 



-n* 
o 



with A = W3 — lj\ — Up — ll>2 — lo\ — u p the one-photon 
detuning, u> p (u> c ) the probe (coupling) laser frequency, 
Clp, c = £ P ,cfJ-p,c/fi the Rabi frequencies associated to the 
corresponding field amplitudes £ PjC and to the dipole mo- 
ments fip^ c of the corresponding atomic transitions. We 
consider an exact two-photon resonance. 

The propagation of two laser pulses in the medium is 
described by the Maxwell equations which, in the slowly- 
varying-amplitude approximation and in the moving co- 
ordinate system 



n 



t — x/c 



(1) 



(such that the original wave operator d/dt + cd/dx be- 
comes cd/drj), read [l| 



drj 



dCl c 
dr] 



(2) 



Here 



27TCiA 



r .j\Hp.c\ 2 N/hc are the coupling factors 
with N the density of atoms. The density matrix el- 
ements pij are determined by the system of equations 
(where the dot denotes the derivative d/dr, since d/dt = 
d /Or) 

pu = HP33 ~ 2Im(0*p 3 i) (3a) 
P22 = 72P33 - 2Im(f2*p 32 ) (3b) 
P33 = -(7i + 12)P33 + 2Im(0*p 31 + Cl* c p 32 ) (3c) 
P31 = -(r + iA)p 3 i + iflp(pn - P33) + iCt c P2i (3d) 
P32 = -(r + iA)p 32 + iCl c {p22 - P33) + iCl p P*>i (3e) 
P21 = ~7cP2i - iCl p pl 2 + i^*P3i (3f) 



III. ADIABATIC PROPAGATION 

It is well known that the use of a counterintuitive se- 
quence of pulses in two-photon resonance, here the cou- 
pling pulse switched on before the probe pulse, leads, in 
the limit of slow evolution of the pulses, to an adiabatic 
passage along the dark state (of null eigenvalue) immune 
to decoherence. The system of linear equations ([3]) can 
be recast as 4-p(es) = L(es)p(es) with L being the asso- 
ciated linear operator and e the formal (small) adiabatic 
parameter (we have omitted here the dependence on rj 
for simplification). The change of variable r = es leads 
to 



e ^( T ) = l ( t )p( t ) 



(4) 



The components of the eigenvector p^ (r) of null eigen- 
value satisfying = L(t)p^°\t) read 



— cos (7, P22 = sin 
p$ = -e^"^ sin I? cos ( 



(0) 

P33 



JO) _ JO) _ 

P31 



(5a) 
4? = (5b) 



with 9 = 9(rj, t) as 



tan/ 



\Qp 



(6) 



In the adiabatic limit, there is thus no coherence between 
the ground states and the excited state, which induce no 
pulse distortion during the propagation [l[ . 

To determine non- adiabatic corrections, we make 
around the adiabatic solution p[r) = p(°> (r) the supera- 
diabatic expansion: 



(o) 
Po 



2 (2) 

e Po 



(J) 



Inserting this expansion in Eq. ([¥]) and identifying the 
orders of e, we obtain a system of linear equations that 
can be solved iteratively. For the first order 



.(0) 

Po 



L Po i 



(8) 



denoting ip — tp p — ip c , we get 



with 71 and 72 the rates of spontaneous emission from 
state 3 to states 1 and 2 respectively, j c the dephasing 
rate between the two ground states, and T the transverse 
relaxation rate. We moreover assume pulses of durations 
much shorter than the characteristic time of decoherence 
of the ground states (usually microseconds in solids) such 
that 7 C will be neglected. The envelopes of the fields 
£ p , c = £p tC (r),T), the Rabi frequencies Cl PtC = Q PtC (r],T), 
and the density matrix components pij = p^ (77, r) de- 
pend on the position 77 = x G [0, L] in the medium of 
length L and on the time r = t — x/c. The bound- 
ary conditions for the Maxwell equations are the pulses 
£ p (0,t) = £ p ,o(t), £ c (0,t) = £ c ,o(t) given at the en- 
trance 77 = x — of the medium. We assume that all 
atoms in the medium are in state 1 before the interac- 
tion with the laser fields. 
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(9b) 

(9c) 
(9d) 

(9e) 
(10) 
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Equations (|9aj) and (|9c|) show that at the first supera- 
diabatic order the excited state stays unpopulated and 
that the adiabatic dynamics along the dark state is sub- 
jected to a loss of order T/il 2 T and to a pure dephasing 
loss of order A/£1 2 T with T the characteristic time of 
interaction. Equations (|9dj) and (|9e|) show the transient 
appearance of coherences between the ground states and 
the excited state, which depend on the first-order non- 
adiabatic loss 9 (and also on ip), but not on the relax- 
ation loss. These coherences will be responsible of pulse 
distortion during the propagation [through the Maxwell 
equations ^] , which will be thus independent of the re- 
laxation loss. Inspection of the solution ([9]) allows one 
to identify the small parameter e <C 1 that permits the 
expansion ([7|) (at least at the first order). We get the 
conditions: 



with 



n 2 T Vt 2 T fl 2 T 
QT»1, - T ;-»l, »1. 

a r 7 C 



(11) 



The two first conditions (large pulse area and small one- 
photon detuning) are the standard condition for adia- 
batic passage for a non-lossy system. The third con- 
dition can be interpreted as the requirement of a long 
interaction time and a large pulse area with respect to 
the decay of the dynamics along the dark state. This 
third condition shows that a fast relaxation FT ^> 1 usu- 
ally encounters in solids can be compensated by a strong 
pulse area satisfying 



(nr) 2 > rr. 



(12) 



We have added in Eqs. (fTTj) the well known fourth con- 
dition considering an additional decohcrence as a pure 
dephasing of rate 7 C between the two ground states. 

It is frequently mentioned that the first condition of 
adiabaticity is not required and that the contr ol p ulse 
may have a rectangular shape (see for instance [11|). It 
is indeed shown in this Ref. [ll[ that in case of a rectan- 
gular control pulse an efficient storage of information is 
possible. However, the parameter T of the above condi- 
tions characterizes the time of interaction when the probe 
and the control pulses both interact, i.e. when they over- 
lap. One can thus consider a rectangular control pulse 
but turned on earlier and turned off later than the probe 
pulse if one preserves the adiabaticity during the field 
overlapping. 

Using the density matrix at first order ([7|) with (J5]) and 
©, we can rewrite the Maxwell equations as 



56> 1 89 
drj u dr 

d_ /n* 

drj \Q 
dip 1 dip 
dij u dr 



(13a) 
(13b) 
(13c) 



Q = Q(r),r) = 
u = u(rj, t) = 



q_ P q c 



q c sin 2 9(i],t) +q p cos 2 9(t],t) 
Q(v,t)S1 2 (t,,t) 



,(14a) 
(14b) 



respectively the two-photon transition strength and the 
group velocity, and the boundary condition 



9(0,t) = 9 (t). 
Denoting the photon fluxes (J — p, c) 



rij = nj(r),T) 



c^t)] 2 _ N\^{ V ,t)\ 



<ij 



i(r),r) = n p (r),r) + n c (r),T), 



(15) 

, (16a) 
(16b) 



Eq. (fl3bl) can be rewritten as dn/dr] = and interpreted 
as the conservation of the photon fluxes during the prop- 
agation: 

„(,,.) = „(0, T , s ^ Ms A r (B^ + &2(l)P). 

V q p Or J 



(17) 

Note that this holds only in the limit of no population in 
the excited state (Le. at the first order of the superadia- 
batic expansion) [10(. 

The obtained system of equations coincides with that 
studied in detail in literature, where pulses of durations 
much shorter than all relaxation times were considered 
(see e.g. [l(| ). However, the system of equations ob- 
tained in previous papers (e.g. in [Ta ]) is valid under 
condition of interaction adiabaticity flT ^> 1, whereas at 
high relaxation rates the inequality (fT2"|) should be ful- 
filled. Solutions of the system of equations (fT3|) have 
been obtained and studied for different regimes of prop- 
agation in a number of works (see, e.g., the review [l[). 
The method of characteristics, recalled in Appendix 
allows one to determine them. 

In particular, equal transition strengths q p = q c lead 
to a constant value Q = q p , and from Eq. (|13b[) to 
f2 independent of 77, thus of shape invariancc during 
the propagation: £1(j],t) = Clo(t — x/c), as well as u: 
u(r),T) = Uo(t — x/c) = £l 2 {t — x/c)/q p with Cto(t) given 
at the entrance of the medium. The pulses [through 
their Rabi frequencies fl p = e ziPp flo(t — x/c) sin 9 and 
il c = e l¥,c f2o(i — x/c) cos 9] propagate in a complicated 
manner through Qq and the mixing angle of solution [see 
Eq. (IA41) ] 



6(v,t)=6(0,£) = 6 (O 



(18) 



where £ = £(x, t) is defined by (|A7[) (where t should be 
replaced by r = t — x/c) [121 ] 



t— x I c 



drfi 2 (r)/g p = x. 



(19) 
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When Slo is additionally constant, the excitation propa- 
gates at constant velocity satisfying 1/u = 1/c + q„/Qn- 
This propagation regime is known as adiabaton [12| . 

In the more general case of unequal transition 
strengths, Eq. (|13a[) entails that the mixing angle still 
reads [13] 



0fa,r)= 0(0,0 = O (£) 



(20) 



and from (THaj) that Q(t),t) = Q(0,0 = Q (0 with 
u(tj,t) = QUO n o(T)/Nq p q c , and £ ee £(x,t) deter- 
mined from Eq. (|A7[) (where t should be replaced by 
t = t — x/c): 



t — x I c 



dr uo(t) 



Q§(0 



(21) 



f2 2 propagates in a more complicated manner than in case 
of equal transition strength as S1 2 («,t) = Qo{O n o{t ~ 
x/c)/N. Since Qo and no do not propagate at the same 
velocity, £7(7?, r) is not any more shape-invariant during 
the propagation. 

If we assume constant phases ip PtC at the entrance of 
the medium, one can conclude from Eq. (|13cP that they 
will thus not change during the propagation. Note that 
this statement is only valid at the exact two-photon reso- 
nance. Evolution of the phase self-modulation for a two- 
photon detuning different from zero has been analyzed in 

In the linear approximation (with respect to the probe 
field), i.e., in the first order with respect to 6, one has 
Q = q Cl fl 2 = f2 2 , u = Sl 2 /q p - We remark that in case 
of equal strengths of adjacent transitions q p = q c , similar 
expressions are obtained with fl 2 in place of fl 2 . Thus, in 
case of equal transition strengths the nonlinearity in the 
probe field results in replacement of the Rabi frequency 
of the control field by the generalized Rabi frequency £l 2 . 
The nonlinearity is more important when the adjacent 
transition strengths are unequal and can lead to forma- 
tion of shock- wave fronts and violation of the interaction 
adiabaticity, as shown in [Toj . 



IV. EFFICIENT STORAGE 



which necessitates a long enough medium. It is usually 
achieved for a strong control pulse switched on before 
and switched off after a weak probe. We present below a 
novel pulse configuration with a strong probe and a weak 
control that allows one to shorten significantly the length 
of storage. The storage requires the stopping of the ex- 
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FIG. 2: Time envelopes of the control (dashed lines) and the 
probe (full lines) Gaussian pulses Qj — exp[— (t/Tj) 2 ], 

j = p, c, at the entrance of the medium with T c — 10T P and 
TT P — 200: a) Usual scheme (shown here for f2 p , ma x/T = 



2.24 x 10 
with Q 



ii. max 



max /r = 2.24 x 10" 

2.24 X 10 - \ fic.ma; 



); b) Proposed scheme 
/r = 2.24 x 10" 2 . 



The storage of the probe pulse is achieved when its 
time profile is mapped into the spatial distribution of 
the atoms in the medium. More precisely, the probe 
pulse can be encoded in the spatial distribution of the 
angle 9o(£(x,t — > +oo)) after the passage of the pulses 
through the medium, inducing for instance a storage in 
the coherence between the lower states (in the adiabatic 
limit): 

P2i(x, t -» +oo) « sin[20 o (£(z, t -> +oo))]. 

(22) 

An efficient storage requires that the excitation, through 
9{x, t), propagates adiabatically until it eventually stops, 



citation, u(x,t) — 0, for all the characteristics connected 
with the initial excitation at the entrance of the medium 
x = 0. We can estimate the length x max required to com- 
pletely store the excitation from (f2"Tj) using the longest 
characteristic, i.e. £ — > — oo and t ~> oo [l(J: 

1 f + °° q 
x m ax = — / dt (g|ft P ,o| 2 + lO^ol 2 ) , 9 = — (23a) 

~ ^ (qnl, max T p + nl max T c ) (23b) 

qp 

with T p and T c the duration of respectively the probe 
and the control, and fij ima x = m a x t^j,o(^)i j — Pi c - 
Minimizing this quantity requires thus q < 1 and a weak 
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probe pulse, as is well known. On the other hand, adia- 
baticity of the interaction is required. Criterion for this 
adiabaticity has been analyzed in , where it has been 
shown that, under the certain condition adiabaticity is 
broken with formation of shock-wave fronts for distances 
exceeding the critical length xq estimated as 

-L p p.max 1 c.max 

x Q v - — '■ . (24) 

<1 P 1-9 

The latter approximation holds for a probe pulse of du- 
ration T p shorter than the control. This shows that re- 
ducing q to lower a; max will however also shorten xq. Pre- 
serving adiabaticity usually requires thus a strong control 
pulse and an efficient storage will take place when 

Zmax < X . (25) 

Through their asymmetry in f2p imax and f^c.max with re- 
spect to q, Eqs. and (f2~4")l show that we can shorten 
x max with the use of a weak control pulse in addition to 
using a small q parameter, while preserving the condition 
of adiabaticity using a strong probe (see Fig. 2 for the 
pulse scheme). 

We have performed numerical calculations of the full set 
of the density-matrix and Maxwell equations for both 
pulses taking into account a fast relaxation TT p = 200. 
It indeed shows that the storage of the probe field in inho- 
mogeneously broadened media is more efficient when we 
use atoms of different adjacent transition strengths. Fig. 
3b presents the temporal evolution of the probe pulse at 
different propagation lengths (normalized as z = xq p /T) 
for the novel scheme of the pulses at the entrance of the 
medium shown in Fig. 2b. The chosen parameters lead 
to f2 c , max T c = £lp,mnxT p w 44, fie )max T c /r « 1, and 
fip max T p /r » 10, which satisfy the adiabatic conditions 
(fn} and (TT2")) at the entrance of the medium (here the 
time of interaction is T ~ T p ). For the comparison we 
present in Fig. 3a the dynamics for the usual scheme of 
pulses used for the optical information storage. One can 
see from Figs. 3a and 3b that in the case of the novel 
scheme the length of information storage is dramatically 
reduced as compared with that of the usual case. Namely, 
for the novel scheme at the propagation length z of about 
1 the probe pulse is completely absorbed by the medium 
while at the same length in the usual scheme the probe 
pulse remains still unabsorbed. 




FIG. 3: Propagation of the probe pulse for different propa- 
gation lengths z for A = 0: a) Usual scheme of pulses with 
Gaussian envelopes at the entrance of the medium displayed 
in Fig. 2a with q — 0.1; b) Novel scheme of pulses (displayed 
in Fig. 2b). 



to reduce the unavoidable losses caused by decoherence. 
The required A systems with very different strengths of 
adjacent transitions are expected to be found in, e.g., 
rare-earth doped solid-state materials. Study to elabo- 
rate efficient schemes to retrieve the stored pulse in this 
strong regime is in progress. 



APPENDIX A: CHARACTERISTICS METHOD 



V. CONCLUSION 

We have shown that the optical information storage 
of a probe in inhomogeneously broadened media of A 
atoms occurs for short medium length when strong probe 
pulses, weak control pulses and very different strengths 
of adjacent transitions are used. For such a scheme the 
optical length may be of the order of unity in contrast to 
the conventional scheme where the optical length must 
well exceed unity. The proposed scheme would allow one 



We briefly recall here the characteristics method (see 
for instance [HI]) to solve the Maxwell equation 

^ + -«0, e(0,t) = 9 (t) (Al) 

in the case of a non-linear group velocity of the form: 
u = u(x,t). We want to transform this linear first order 
partial differential equation into an ordinary differential 
equation along the characteristic curves (x(s),t(s)). Us- 
ing the chain rule £e(x(s),t(s)) = f§|f + § f s , we ob- 
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tain 



d_ 

ds 



0(x(s),t(s)) = 0, 



(A2) 



if we set (choosing x(0) = and denoting £ := t(0)) 
dx 



1, i.e. x(s) = s 
ds 

dt 1 . 
— = -, i.e. t(s) 



ds' 



(A3a) 
t (A3b) 



ds u J u(x(s'),t(s')) 

Equation (|A2[) leads to the solution 

8(x(s), t(s)) = 8(x(0), t(0)) - 0(0, - O (£) (A4) 

which is thus a constant along the characteristics. £ al- 
lows one to label a characteristics i(x) rewritten from Eq. 
(|A3b|) as (since s = x) 



t(x) 



/ u(x',i(x')) 
The non- linear velocity u derived from (jA5|) : 



1 _ rft 
u dx 



(A5) 



(A6) 



corresponds to the instantaneous velocity along a char- 
acteristic. 

When u is independent of t, i.e. it = u(x), one has 
£ = t — J Q dx' /u(x'), which becomes £ = t — x/u if it is 
also independent of x. 

When u depends only on t, i.e. u = u(t), Eq. (|A3b|) 
can be rewritten as 



o?i' it(i ) = s = x, 



which gives in principle £ = £(x, i). 



(A7) 
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